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We investigate quark deconfinement by calculating the effective potential of the Polyakov loop 
using the non-perturbative propagators in the Landau gauge measured in the finite-temperature 
lattice simulation. With the leading term in the 2-particle-irreducible formalism the resultant ef- 
fective potential exhibits a first-order phase transitions for the pure SU(3) Yang-Mills theory at 
the critical temperature consistent with the empirical value. We also estimate the thermodynamic 
quantities to confirm qualitative agreement with the lattice data near the critical temperature. We 
then apply our effective potential to the chiral model study and calculate the order parameters and 
the thermodynamic quantities. Unlike the case in the pure Yang-Mills theory the thermodynamic 
quantities are sensitive to the temperature dependence of the non-perturbative propagators, while 
the behavior of the order parameters is less sensitive, which implies the importance of the precise 
determination of the temperature-dependent propagators. 

PACS numbers: 11.30.Rd, 12.40.-y, 21.65.Qr, 25.75.Nq 



Introduction It has been a long-standing question in 
physics of quantum chromodynamics (QCD) how to un- 
derstand confinement of quarks and gluons in the vac- 
uum and how to clarify the nature of deconfinement in 
a medium in extreme environments (see Ref. [1] for re- 
views) . 

It was Polyakov [2] who first addressed the deconfine- 
ment phase transition successfully in the strong-coupling 
limit of a pure Yang-Mills theory. The order parameter 
for deconfinement was then identified, which is now called 
the Polyakov loop. Later on, the strong-coupling expan- 
sion was extended to implement quarks and the chiral 
dynamics [3]. One of the most popular approaches in 
the QCD phase-diagram research, i.e. the chiral effective 
model such as the Nambu-Jona-Lasinio model [4-7] , the 
linear-sigma model [8], etc [9] with the Polyakov loop, is 
a natural extension from the strong-coupling QCD. 

The largest ambiguity in the P-chiral models lies in 
the choice of the effective potential of the Polyakov loop. 
Though the initial choice was motivated by the strong- 
coupling expansion [4], it is now common to adopt the 
potential that fits the pure Yang-Mills thermodynamics 
from the lattice simulation either in a polynomial form [5] 
or in a Haar- measure form [4, 6]. Since the fitting pro- 
cedure does not refer to microscopic dynamics at all, it 
is unclear how the Polyakov-loop potential is related or 
unrelated to non-perturbative characteristics near T c as 
studied in the matrix model [10] (see also Ref. [11] for an 
insight into the physical contents of the potential). 

An important breakthrough came from an attempt 
to understand quark deconfinement in terms of the 
Landau-gauge propagators that describe gluon confine- 
ment [12, 13]. In the Landau gauge it is the deep-infrared 
enhancement in the ghost propagator that causes confine- 
ment, while the gluon propagator is infrared suppressed. 
This behavior is qualitatively consistent with the confine- 
ment scenarios by Kugo and Ojima (they are equivalent 



if the ghost renormalization factor diverges at zero mo- 
mentum [14]) and also by Gribov and Zwanziger [15]. 
Indeed the gluon and ghost propagators in the Landau 
gauge at zero and finite temperature have been calculated 
in the lattice simulation [16-18], the Dyson- Schwinger 
equation (DSE) [19, 20], the functional renormalization 
group (FRG) approach [13, 20, 21]. 

In this work we report on an extension of Ref. [13] 
using the T-dependent propagators from the lattice- 
QCD simulation [18] (see Ref. [22] for the finitc-T 
propagators from the FRG calculation). Also, we try 
to address the thermodynamic quantities constructed 
from the propagators. We would stress here that our 
central point is not only an extension of Ref. [13] but 
it should aim to build a bridge over the model studies 
and the first-principle functional approaches. Such a 
work must be extremely useful for both sides; there 
are many arguments to suggest that the back-reaction 
from the matter to the glue sector has crucial impacts 
on the QCD phase-diagram research based on the 
P-chiral models [8, 23]. In principle this shortcoming of 
insufficient back-reaction in the P-chiral models could 
be resolved with the DSE or FRG approaches [24], but 
then one has to confront huge-scale computations. As 
we see later, thanks to a nice parametrization of the 
propagators in Ref. [18], our practical strategy requires 
only minimal modifications in the model-study pro- 
cedures. The outcomes are quite promising as we will see. 

Propagator and Parametrization In the covariant 
gauge the gluon propagator inverse can be decomposed 
into the transverse and the longitudinal parts using the 
projection operators; 



P^Pv 



(1) 



2 



T 


r? (GeV 2 ) 


dt (GeV -2 ) 


rf (GeV 2 ) 


di (GeV~ 2 ) 


0.86T C 


0.880 


0.143 


0.256 


0.220 


1.20T C 


0.963 


0.140 


1.018 


0.162 



TABLE I. Parameters for the non-perturbative propagators 
at T = 0.86T C and 1.20T C taken from Ref. [18]. 



and the further decomposition relative to the rest frame 
at finite temperature, 

Pi, = (i-W(i-M(v-^). p£v = T^-P*. 

(2) 

The gluon propagator in the Landau gauge is then 
parametrized as 

Da^P 2 ) = 5 ab (D ( PP^ + D ( T L) P^ + ZD L L^) . (3) 



Interestingly these propagators as well as the ghost prop- 
agator Dc are compactly expressed in a Gribov-Stingl 
form; 



D 
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dt( P 2 + d^) 

(p2 + r 2)2 
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Dl = — , 

p 2 

D C = 



A" 1 



(4) 



(p2) 2 



as discussed in Refs. [18, 25]. This form also ap- 
pears in a low-energy effective description of the 
Yang-Mills theory [26]. We note that we postulated 
non-renormalization for the longitudinal gluon propaga- 
tor because Zl = 1 + C(£) and £ — > is taken after all 
in the Landau gauge. In this work we use two sets of 
parameters at T = 0.86T C and T = 1.20T C taken from 
Ref. [18], which is listed in Tab. I. The ghost propagator 
can be well reproduced by a choice of d~ 1 = 0.454 GeV 2 
which is nearly independent of T. 

Polyakov Loop Potential The Gribov-Stingl form is 
especially convenient for the practical computation of 
the partition function and the effective potential of the 
Polyakov loop. In the present study we make use of an 
approximation motivated from the 2-particle-irreducible 
(2PI) formalism or the FRG equation. In the 2PI formal- 
ism for example, the full effective action can be expressed 
generally as 

r = itrlnG- 1 - itr ln(G _1 -G 1 )G + T 2 [G] (5) 

with the full propagator G, the tree- level propagator Go, 
and the sum of the 2PI diagrams T 2 [G] . Thus, all known 
approximations are derived from some truncation onto 
Eq. (5); the Hartree approximation for example picks up 
the first term trlnG -1 only with an assumption that G 
is written with a mean-field mass. Usually the full propa- 
gator should be fixed self-consistently within the 2PI for- 
malism. Now that we already know the non-perturbative 



propagators, however, we can directly plug them into 
Eq. (5) to have better results. In this "hybrid" method, 
the dominant contribution to the gluonic pressure comes 
from 



/3Q g lue - -itrlnZ^ 1 +trlnD c ; 1 



= --trln£>; 



-trln D 



(T)-l 



- -tr In D 



(L)-l 



+ trlnL> 



c 



(6) 



We note that this leading-order truncation corresponds 
to the so-called quasi-particle model of thermodynam- 
ics [27]. In other words this method presumably works 
in the deconfinement regime where gluons are the physi- 
cal degrees of freedom, while the approximation may fail 
the regime of glucballs at low T. Importantly, in prin- 
ciple, one can systematically improve the approximation 
by evaluating the sub-leading terms and T 2 diagrammat- 
ically once Da and Dc are known. 

For the purpose of calculating the effective potential for 
the Polyakov loop, we keep a constant A4 in the temporal 
covariant derivative; p 2 in the gluon and ghost propaga- 
tors are replaced with p 2 — {2-kTti + gfiAi) 2 + p 2 in the 
color adjoint representation. We note that the Polyakov 
loop is defined as 



^trVexp^ig J A 4 dx^j , (7) 



$ = -trL 3 = 
3 3 3 



where V represents the time ordering and A4 is a matrix 
in the color fundamental representation. 

Then, because the propagators (4) are written as a 
combination of (p 2 +m 2 ) -1 , it is straightforward to carry 
out the summation over the Matsubara frequency in an 
analytical way, i.e. for the transverse gluons for example, 
wc have 



trlnD 



(T)-l 



2trln(p 



z2 



- trln(p 2 + d t 1 ) 



= 2W B {r 2 t )-W B {dt 1 ) 



(8) 



where we have defined, 
W B (m 2 ) = -2V J -^trln(l-L 8 e-^ / P 2 +™ 2 ) (9) 

with L s being the Polyakov loop matrix in the color 
adjoint representation; (Ls)a6 = 2tr(i a L 3 i;>Lg). We 
dropped the divergent zero-point energy that is indepen- 
dent of T and thus an irrelevant offset. One can take the 
trace in color space explicitly to find [11], 

W B (m 2 ) = -2V J ^ln(l + £c„e-W^) , 

(10) 

where G 8 = 1, C x = C 7 = 1 - 9$$, G 2 = G 6 = 
1 - 27$$ + 27($ 3 + $ 3 ), G 3 = G 5 = -2 + 27$$ - 

81($$) 2 , G 4 = 2 [-1 + 9$$ - 27($ 3 + $ 3 ) + 81($$) 2 . 
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FIG. 1. Polyakov loop potential near T c that exhibits a first- 
order phase transition as expected for the SU(3) case. With 
the parameter sets at T = 0.86T C and T = 1 .20T C in Tab. I 
we find T c = 289 MeV (as shown) and T c = 351 MeV, respec- 
tively. 



In this way, using Eqs. (6), (8), and (10), we can nu- 
merically calculate the Polyakov-loop effective potential 
/30gi ue [$] , and then determine $ from 9fi g i ue [$]/9$ = 0. 
We find that the critical temperature in our treatment is 
T c = 289 MeV and T c = 351 MeV for the parameter sets 
at T = 0.86T C and T = 1.20T C in Tab. I, respectively. 
If we utilize the full T-dependent propagator, thus, the 
critical temperature should be somewhere in between. 
Considering the empirical value T c ~ 280 MeV [28], one 
might wonder that our estimate of T c is a bit too larger. 
This is fine because neglected 1^ has an effect to push 
T c down [29]. This semi-quantitative agreement of T c is 
amazing for its simplicity. The energy scales are provided 
through r 2 and d as given in Tab. I and there is no ad- 
justable degrees of freedom. It is also worth mentioning 
that the Polyakov loop potential formulated here leads to 
a second-order phase transition for the color SU(2) case. 

Although T c is such different depending on the param- 
eter sets, interestingly, $ and all thermodynamic quanti- 
ties as functions of T are identical if plotted in the unit 
of T c . This is a quite non-trivial finding, and supports 
the validity of the fitted Polyakov loop potential charac- 
terized by only one dimensionful parameter To [5, 6]. 

We show thermodynamic quantities as a function of 
T in Fig. 2 where the normalized pressure and the in- 
teraction measure are plotted. Because we adopt the 
parameter sets at T = 0.86T C and T = 1.20T C (yielding 
indistinguishable results if scaled with T c ), we choose the 
temperature range up to 1.3T C here. We make a compar- 
ison with the lattice data and confirm semi-quantitative 
agreement. In particular, as seen in Fig. 2, the agree- 
ment of the pressure is pretty good, while the interaction 
measure does not fit the lattice data well. 

The discrepancy above ~ 1.2T C should be attributed 
to the neglected T-dependence in the non-perturbative 
propagators as seen in Refs. [18, 31]. Once it is correctly 
taken into account, we expect that the overshoot of 




FIG. 2. Thermodynamic quantities as compared to the lattice 
results taken from a table in Ref. [30]. The numerical results 
are indistinguishable for the parameter sets of T = 0.867^ and 
T = 1.2T C . 



thermodynamic quantities could become milder. It 
should be mentioned that, in the T oo limit, ther- 
modynamic quantities in our method approach the 
Stefan-Boltzmann limit, by construction, as all propa- 
gators go to ~ 1/p 2 . Below T c , on the other hand, we 
find that some of thermodynamic quantities go negative, 
which is caused by too strong ghost contributions. It 
is an unanswered question how to extract the expected 
behavior of the glueball gas [11, 30] from the gluon 
and the ghost propagators. One should cope with 
T'2 in Eq. (5) in the confinement regime. Although 
one may think that glueballs are too heavy to make 
a sizable contribution to thermodynamics, the electric 
glueballs can be significantly light in the vicinity of the 
(second-order) critical point [32]. 

Dynamical Quarks Just for the demonstration pur- 
pose of the usefulness of our effective potential let us 
apply it to an effective model. We emphasize that our 
goal is not to analyze the model itself but to seek for 
the possibility to utilize this inverted Weiss potential 
in phenomenology that is complementary to the field- 
theoretical argument [33]. To this end we adopt the co- 
variant coupling in the (2+l)-flavor quasi-quark descrip- 
tion. Then, the thermodynamic potential from the quark 
contribution reads, 



/^^quark 



6 ^E f^p 2 + M f-^ w ^ M 



+ g s ((uu) 2 + (dd) 2 + (ss) 2 ) + 4g d (uu) (dd) (i 



(11) 



where M u = m u — 2g s (uu) — 2gd{dd) (dd) etc and the last 
two terms above represent the condensation energy. Here 
we defined, 



W F (M 2 ) = V 



d 3 p 
(27p 



trln 1 



(12) 
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FIG. 3. Order parameters (see Eqs. (13) and (14) for defini- 
tion) from our model calculation and the lattice-QCD results. 
The bands in our results correspond to the uncertainty from 
the parameter choices at T = 0.86T C (with T c = 182.5 MeV 
for the unit of the horizontal axis) and T = 1.20T C (with 
T c = 191 MeV). The gray bands are the lattice-QCD re- 
sults in the continuum extrapolation from Refs. [35] (with 
T c = 156 MeV). 



in the color fundamental representation. Thus, we 
choose the parameters according to the standard set 
of the (2+l)-flavor NJL model as A = 631.5 MeV, 
g s A 2 = 3.67, g d A 5 = -9.29, m„ = m d = 5.5 MeV, 
m s = 135.7 MeV [34]. 

We can express Wf(M 2 ) in terms of $ [5] and solve the 
Polyakov loop <£> and the chiral condensates (qq) to min- 
imize the total potential Q g \ ue + f2 qua rk- Figure 3 shows 
our numerical results for the order parameters; The nor- 
malized chiral condensate is defined as 



A 



(uu) - (m ud /m s )(ss) 
(uu) - {m ud /m s )(ss) 



(13) 



where (uu) = (dd) represents the light-quark chiral con- 
densate and the denominator is the value at T = 0. We 
adopt the strange-quark number susceptibility, 



/T 2 = J_d^ 
Xs/ T 2 dfi s 



i o 2 n 
vt 2 ai4 



(14) 



for the deconfinement order parameter instead of the con- 
ventional Polyakov loop. This is because the Polyakov 
loop has a large renormalization factor, and in the mean- 
field approximation $ should be considered as an in- 
ternal parameter rather than an observable. To draw 
Fig. 3 we chose T c = 182.5 MeV for the parameter set 
at T = 0.86T C , T c = 191 MeV for that at T = 1.20T C , 
and the lattice-QCD results from Ref. [35] are plotted 
as a function of T/T c with T c = 156 MeV. One might 
have an impression that T c in the model side is too large 
as compared to the lattice-QCD data, but this is to be 
improved with the back-reaction from the polarization 
diagrams [8, 23]. 

In this case with dynamical quarks thermodynamic 
quantities do not show unphysical behavior near T c be- 
cause quark degrees of freedom dominate over gluons. 
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FIG. 4. Thermodynamic quantities as compared to the 
lattice-QCD results taken from Ref. [36]. The bands cor- 
respond to the uncertainty from the parameter choices at 
T = 0.86T C and T = 1.20T C . The choice of T c is the same 
as in Fig. 3. 



The pressure and the interaction measure are plotted in 
Fig. 4. We see that our numerical results quantitatively 
agree with the lattice-QCD data taken from Ref. [36]. 

Our method with the momentum integration in 
Eq. (9) is simple enough to be an alternative of the 
Polyakov-loop potential used in the market of the 
P-chiral models. Moreover, it is advantageous for 
our method to be extendable to implement missed 
contributions such as the screening effects through the 
quark polarization diagrams. The improvement in this 
direction should be extremely important to tackle the 
realistic QCD phase structure especially at high baryon 
density or strong magnetic field [37] (for a recent review, 
see Ref. [38]). These effects do not directly couple to 
gluons, and nevertheless, the nature of deconfinement is 
affected substantially through the quark loops that carry 
the baryon number and the electric charge. Progresses 
in this direction shall be reported elsewhere. 

Summary We elucidated how to construct the ef- 
fective potential of the Polyakov loop from the non- 
perturbative propagators of gluon and ghost in the Lan- 
dau gauge available from the lattice simulation. This 
is an extension of the idea of Ref. [13]. We took the 
fitting forms of the finite-temperature propagators from 
Ref. [18] and found a quite tractable way to calculate 
thermodynamic quantities as well as the order parame- 
ters as functions of T. We showed that the thermody- 
namic properties are nicely consistent with the lattice 
data in the vicinity of T c . Furthermore we introduced 
dynamical quarks in the quasi-quark approximation to 
reproduce the simultaneous crossovers of deconfinement 
and chiral restoration. We made sure that our poten- 
tial works well even on the quantitative level without 
fine-tuning of any parameter. It would be an intriguing 
future problem to apply our Polyakov-loop potential to 
the non-local version of the chiral model [391. 



■5 



We believe that this present work takes one step for- 
ward to the understanding of the QCD phase diagram in 
extreme environments based on the first-principle-type 
calculations. 
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